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ABSTRACT 

We investigate the production of the lightest neutral Higgs boson in as- 
sociation with tau sleptons in the CP violating minimal supersymmetric 
standard model (MSSM) at future high-energy e + e~ linear colliders. In pa- 
rameter space of the constrained MSSM, the production cross section of 
e + e~ — > h°TiTi can be very substantial at high energies. This process 
would provide a new production mechanism of neutral Higgs bosons and 
open a window to probe the Higgs-stau coupling and some soft supersym- 
metric breaking parameters at next linear colliders (NLC). 
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I. INTRODUCTION 



Supersymmetric (SUSY) theories, in particular the MSSM [[T], ||, are cur- 
rently considered as the most theoretically well motivated extensions of the 
standard model (SM). It provides a natural solution to the hierarchy prob- 
lem ||, the generation of the electroweak symmetry breaking @], the grand 
unification of the gauge coupling constants ||, as well as the elucidation 
of the baryon asymmetry and the cold dark matter in the universe || |7j. 
The MSSM contains five physical Higgs bosons ||: two neutral CP-even 
scalars (h° and H°), one CP-odd neutral scalar A , and a pair of charged 
scalars {H + , H~). Such particles break the electroweak symmetry and are 
responsible for the masses of the fermions and bosons. The search for the 
Higgs bosons and the study of their properties are one of the major goals 
of present and future colliders 0. Contrary to hadron machines, at e + e~ 
colliders the background is relatively clean. A future high-energy and high- 
luminosity e + e~ linear collider operating at energies \fs = 500 ~ 1000 GeV 
would provide an ideal environment for searching the lightest Higgs boson 
and performing precision measurements of it [IIJ . 

In the context of the MSSM, some typical production mechanisms of 
Higgs bosons, such as Higgsstrahlung [II]], vector bosons fusion flI5| , Higgs 



pair production [13|], and associated production of Higgs boson with top or 
stop [|14j], have all been extensively studied at NLC energies. In addition, 
it was realized that couplings of Higgs bosons to sfermions can also provide 
additional sources for Higgs bosons production [|lj|, in particular, when the 
mixing between the third generation scalar taus is very large, the associ- 
ated production Higgs bosons with tau sleptons might provide a important 
information on the soft SUSY breaking scalar potential. In this paper, we 
will discuss in more detail the last type of the prospecting processes in the 
CP violating MSSM. Our motivation is as the following: First, for the large 
Yukawa coupling of the tau slepton, the mixing between tau sleptons can be 
very large. As a result, the lighter stau can not only be significantly lighter 
than the top quark and all scalar quarks, but also at the same time their 
couplings to Higgs bosons can be strongly enhanced. Therefore, the reac- 
tion e + e~ — > h°T^r^ can be more phase-space favored, and also it would 
be accessible to produce and measure them at the NLC. Second, it has re- 
ceived growing attention that CP-violating phases associated with the third 
generation trilinear soft term can induce CP violation in the Higgs and the 
third generation sfermion sectors through quantum corrections ]16j. Such 
phases may allow baryogenesis and do not necessarily violate the stringent 
bound from the nonobservation of electric dipole moments (EDM) |T7|, I8|. 
In fact, some of these phases can be 0(1), so as to provide non-SM sources 
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of CP violation required for dynamical generation of the baryon asymmetry 
of the universe |1| . Since these phases directly affect the couplings of Higgs 



bosons to third generation sfermions, moreover, these couplings play an im- 
portant role in the phenomenology of the third generation sfermions and 
Higgs bosons at colliders, we will consider these effects of the CP-violating 
phases on the associated Higgs production with the stau*-stau pair. Last, for 
most of the SUSY parameter space allowed by present data constrains, the 



MSSM Higgs sector is in the decoupling regime p0[ . In this scenario, only 
the lightest Higgs boson would be accessible at the NLC, final states with 
the other Higgs bosons or with r^Tvf pair should be phase-space suppressed. 
In what follows, therefore, we will focus on the case of the lightest h° bo- 
son of MSSM in the decoupling regime, and discuss only the production in 
association with lighter tau sleptons. 

The remainder of this paper is organized as follows. In Section II we 
outline the masses and mixing of the tau sleptons, and then we present the 
relevant couplings and analytical expressions of the production cross section 
for the process e + e" — ► hpT^r^. In Sec. Ill, a detailed numerical analysis of 
the cross section is given in an appropriate MSSM scenario satisfying the ex- 
perimental and theoretical constraints. Sec. IV is devoted to the conclusions. 

II. MASSES, COUPLINGS AND CROSS SECTION 

In this section, first, to fix notation we will simply summarize the masses 
and mixing of the tau slepton sector in the MSSM, and then list the relevant 
couplings involved in the associated production of the lightest neutral Higgs 
boson with tau sleptons. Finally, we will give analytical results of the cross 
section of e + e~ — > hPr^r^ production. 

A. Tau slepton mass and mixing 

In the tau slepton sector jl], 0, the left- and right-handed tau sleptons 
can strongly mix, the mass-squared matrix is given by 

2= f M 2 L m T (A T + fittm[3) \ 

? \m T (A* T + fi* tan (3) M\ J { ) 

with 

M 2 L = m 2 ~ L + m 2 T + {\m 2 z - M 2 W ) cos 2(3 , 

(2) 

M\ = m% R +m 2 T + {M 2 Z - M 2 W ) cos 2(3 , 
where the left- and right-handed soft SUSY-breaking tau slepton 
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masses, respectively. The soft breaking trilinear stau coupling A T and Higgs 
mass mixing parameter /i are complex, 

A T = \A r \e iv , fi=\fi\e lT1 . (3) 

The phases (p and r\ are the sources of CP violation. In fact, the phase rj can 
be transferred into the phase (p in the form of (p — t] by the field redefinition. 
So we will take /i to be real, and vary (p from to 2tt in the following. The 
tau slepton mass eigenstates can be realized by a unitary transformation U 
which diagonalizes the mass-squared matrix M~ , 

UM*rf = diag(m| l , mfj , (4) 

where the diagonalization matrix can be parameterized as 

tt _ ( cos9 T sm9 T e t5 \ 

\ -sm9 T e- iS cos9 T J ^ 

with 

5 = arg (A T + /itan/3) . (6) 
The tau slepton mixing angle and mass eigenvalues are then given as 
2m T \A T + /itan/3| 

(7) 



tan 26* 



Ml -Ml 



1 

™n, 2 = 2 



Ml + M^T J (M 2 L - M 2 R ) 2 + 4m2|A r + /xtan/?| 2 



with convention ^ ^ 7r/2, m| ^ m~ . For large values of |^4 T | ; l/^l and 
tan/3, the mixing in the tau slepton sector can be very strong. In the con- 
strained MSSM or minimal supergravity (mSUGRA) scenario |J1 [2~H , the soft 
SUSY breaking scalar fermion masses and gaugino masses are respectively 
unified as m and m\ii at the grand unified theory (GUT) scale Mqut- The 
left- and right-handed tau slepton masses at the weak scale 0(Mz) are then 
given in terms of the universal masses and m\ji through running renor- 
malization group equations (RGE) as 



m~. ~m„ + 0.52 m\ /2 - 0.27M| cos 2/3 , 

m~ R ~m2 + 0.15m? /2 -0.23M|cos2/3. ( ^ 

Inputting values of the parameters \A T \, (p, tan/3, m , mi/ 2 , the mass 
spectra and mixing angle of the tau sleptons can then be worked out. For 
smaller values of m and mi/2 , the stau masses are lower. 
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FIG. 1. Feynman diagrams for e + e — > h G T^r l production. 

B. e + e — > h°Tj" r{" production 

e + e~ — > h Ti~Ti production proceeds via three kinds of diagrams, as 
shown in Fig. 1. When the h°T*Ti couplings are large, the most important 
is Higgs boson emission from the t\ states which are produced through s- 
channel photon and Z-boson exchange (Diag. a and b in Fig. 1). One also 
has the conversion of t% to ri, which involves the h°T*T2 vertex (Diag. c and 
d in Fig. 1). The last type is h°Z* production with Z* — > rf (Diag. e in 
Fig. 1). Though the two last types of Feynman diagrams turn out to give 
small contributions, there are some regions in the parameter space where they 
cannot be neglected. The relevant couplings for the process are described by 
the following set of interaction Lagrangian in the MSSM [|l], §]: 
(a). The couplings of a neutral Higgs boson h° to a pair of tau sleptons 
TiTjfaj = 1,2) is given by 

£fc°7?53 = 9ijh Q T*Tj (9) 
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with 



9ij 



gm z 



c w 

gm z 



i. s w ~ o) sin ( a + 0) + 



sina 



m z ) cos (3 



+ 

+ 
+ 



c w 
gm T 

gm T 

2c w m z 



-s w sin(a + (3) + 



A 



sin a 

T cos/5 
sin a 



m z 
cos a 



sm a 



cos /3 
cos a 



cos j3 cos /3 



cos (3 

(4)*4 
(4)*4 , 



(4)H 
(4)*4 



(10) 



where sy^ = sin9w,cw = cos 9w, o> is the CP-even Higgs mixing angle at 
the tree level. When the loop corrections are included, the tree level pa- 
rameters are deformed, even inducing the mixing among Higgs with different 
CP properties p3 |. In the decoupling limit and up to very small radiative 
corrections, however, the mixing angle a reaches the limit (3 — n/2. 
(b). The couplings of neutral vector gauge bosons (7 and Z°) to tau slepton 
pairs are written as 



W_ 

c w 



swewSijA* + (-s^Sij + ^(4)*^/) Z» 



(c). The couplings of neutral vector gauge bosons to charged lepton pairs 
are described by the same one as in the SM 



C 



Vff 



_g_ 

c w 



ls w cwA»fo,f + Z»f llM 



_ o 2 



1 -7 5 2 1 + 7 5 



/ -(12) 



(d). The couplings of a neutral Higgs boson to a pair of massive neutral 
gauge bosons is given by 



C 



h°ZZ 



gm z 

2c 



w 



sin(/3 - ^hPZ^Z^ . 



(13) 



Using Feynman rules, the Dalitz plot density of the production cross section 
can be obtained. In terms of the energies Ex, E 2 , E^ of the three final state 
particles r-f, t+ and h°, respectively, it is written as 



da 



dExdE 2 8(2 



+ 



7T 



771 



s 2 — Asm 2 1 8 



s -Ap + ^^Re(AB') 



HI- ^ + {El - ml) te + 



\B\ 



(14) 
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where as usual s is the total center-of-mass energy squared, the coefficients 
A and B are defined as follows: 



.4 
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cwmz 



\U- 



9n 



t 12 

nl 




+ 



2y/sE 2 



+912 



+921 



{u\ 2 yu : 



s - 2^~sE 2 + m\ 



m% 
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s — 2^fsE\ + m\ — m\ 



B 




' w 



— s 



g sin(/3 - a) (s - y/sE h - m|) 
c w m z {s - m|) 

g 3 sin(/5 — a) 



<? w m z {s - m|) 



(15) 



To obtain the total production cross section er(e + e — > h ^^ ), one can 



rewrite equation (14) in the standard form for the Dalitz plot density |24 
and then integrate over the variables si 2 , S23 whose integration bounds are 
given from kinematics. For a given center of mass energy the cross section 
depends on the masses of created particles and vanishes at the kinematic 
boundary 

III. NUMERICAL RESULTS 

In this section, we will illustrate our numerical results of the total produc- 
tion cross section for e + e~ — > hPr^T^ based on a specific constrained MSSM 
scenario for the relevant SUSY parameters. The whole parameters appearing 
in our analyses contain \A T \, ip, tan/?, m , mi/ 2 , mh- Including the two- 
loop radiative corrections, the mass of the lightest Higgs boson, m^, depends 
mainly on the following variables [ 2~5| : the top quark mass, the top and bot- 
tom squark masses, the mixing in the stop sector, the pseudoscalar Higgs 
mass mA and tan j3. In particular scenarios, like in SUGRA or gauge and 
anomaly mediation |21], |26|, |27| , some of the soft parameters may be related 
to each other. However, since the mechanism of SUSY breaking is unknown, 
the soft SUSY breaking parameters associated with the tau slepton sector 
can be considered to be independent of the parameters in the Higgs and 
squark sector P8|| . in what follows, therefore, we consider the former seven 
parameters as free phenomenological parameters to be constrained by exper- 
imental and theoretical considerations. Among the constraints that we are 
going to impose on the MSSM are those which follow from the comparison of 
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FIG. 2. The cross section cr(e + e~ — > /^r+r-f ) as a function of the scalar 
fermion mass mo for tan/3 = 5, = 500 GeV, |A T | = 1 TeV, tp = 7r/4, 
m/j = 120 GeV. Three curves correspond respectively to three representative 
values of the gaugino mass m 1 / 2 ■ 



the SM with experimental data, from the experimental limits on the masses 
of yet unobserved particles , etc, and also those that follow from the ideas 



of unification and from SUSY GUT models [^9 . 

In the following numerical analysis, we fix the center of mass energy to 
y/s = 1 TeV, and then demonstrate in turn the dependence of the cross 
section on the choice of parameters. In Fig. 2 and Fig. 3, the cross section 
is shown as a function of the unified scalar fermion mass m for the case of 
low and high tan/? value: tan/3 = 5 and tan/3 = 40, respectively. Three 
curves in each figure correspond respectively to three chosen representative 
values of the unified gaugino mass m^j . The other parameters are typically 
taken, for example, as the following: \A T \ = 1 TeV, ip = tt/4, = 120 
GeV, \fi\ = 500 GeV (for tan/3 = 5), = 300 GeV (for tan/5 = 40). As 
can be seen, a typical magnitude of the cross section is 0.002 ~ 0.03 fb. 
The variation of magnitude depends mainly on m . As the value of m is 
increased, the cross section drops rather dramatically. The values of m 1( / 2 
do not change the whole trends of the plots, but can slightly shift the value 
of the cross section in the same order of magnitude. This is due to the fact 
that the tau slepton masses basically depend on m . In addition, the value 
of tan j3 can also influence the value of the cross section obviously, which can 
be seen from the position interchange of the solid line and the long-dashed 
line in the two graphs, in particular, for the large tan /3, smaller mo and 
are excluded by experimental limits. As the value of tan/? increases from 
low to large, the cross section decreases significantly. However, for no more 
than about 200 GeV value of m in the case of the low tan (3 (about 150 GeV 
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FIG. 3. The cross section cr(e + e~ — > hPT^Ti) as a function of the scalar 
fermion mass mo for tan/3 = 40, |/i| = 300 GeV, \A T \ = 1 TeV, ip = 7r/4, 
m^ = 120 GeV. Three curves correspond respectively to three representative 
values of the gaugino mass m 1 / 2 ■ 



in the case of the large tan/3), the production cross section can exceed the 
value a(e + e~ — > /^r+r-f ) ~ 0.01 fb. With a integrated luminosity 500 fb~ , 
this provides several dozen events in a few years, which should be sufficient 
as a sample to isolate the final state and measure the gh°r*Ti coupling with 
some accuracy. 

In the Fig. 4, the cross section is displayed as a function of the Higgs 
mass mixing parameter for three representative values of the lightest 
Higgs boson mass rrih , and the other parameters are fixed to tan /3 — 5, 
| A,- 1 = 1 TeV, tp = 7r/4, m = 200 GeV, m 1/2 = 250 GeV. The value of \fi\ is 
constrained by the requirement of radiative electroweak symmetry breaking 
P|, while the values of m^ are chosen from the current experimental and 



theoretical considerations |25|, [3(J. The plots show that the cross section 
changes somewhat large with values of but is not too sensitive to the 
values of m^. For a larger value of \fi\, the cross section is smaller. When the 
value of rrih is increased by 20 GeV, the cross section drops approximately 
0.001 fb. 

The CP violating phases in A T , arg(A r ), can modify the masses and the 
mixing parameters of stau as well as the Higgs-stau*-stau couplings. This 
would also have large impact on the production cross section of the hPr^r^ 
final state since it is directly proportional to the square of the h^r^Ti cou- 
plings. In Fig. 5, the cross section is shown as a function of the CP-violating 
phase (p (=arg(v4 r )) for three representative values of the soft breaking tri- 
linear stau coupling \A T \ , and the other parameters are fixed to tan/3 = 5, 
\fi\ = 500 GeV, m h = 120 GeV, m = 200 GeV, m 1/2 = 250 GeV. The graph 
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FIG. 5. The cross section <r(e + e~ — > h°T^T^) as a function of the CP- 
violating phase tp for tan/3 = 5, = 500 GeV, rrih = 120 GeV, mo = 
200 GeV, TO1/2 = 250 GeV. Three curves correspond respectively to three 
representative values of the soft breaking trilinear stau coupling \A T \ . 
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shows that the cross section has a symmetry about tp — n. This results from 
that the stau mass are symmetric under tp — > tp — 2n. Furthermore, the 
cross section reaches a maximum near tp — ir. This is because the coupling 
ghPr^n has a maximum at ip = ir and its phase dependence is stronger than 
that of the masses of stau. However, the trends of plots and the varying 
magnitude of the cross section are obviously subject to the values of \A T \, in 
particular, it should be noted that the cross section might increase or also 
decrease with the values of \A T \ in deferent interval of the phase tp. 

IV. CONCLUSIONS 

In summary, in the framework of the CP violating MSSM, we have dis- 
cussed the production of the lightest neutral Higgs particle in association 
with the lighter tau slepton pair, e + e~ — > h°r^r{ , at future high-energy 
e + e~ linear colliders. Simple analytical formulae for the cross section of the 
three-body process have been given. In a specific constrained MSSM sce- 
nario, where the lightest Higgs boson is in the decoupling regime, we have 
analyzed in detail the dependence of the cross section on the relevant SUSY 
parameters. On the one hand, for large values of \A T \, and tan/5, the 
mixing between tau sleptons can be very large, this results in the lighter stau 
can not only be rather light but also at the same time their couplings to Higgs 
bosons can be strongly enhanced. As a result, the cross section can be very 
substantial. On the other hand, the effect of the CP violating phase in A T 
on the cross section is also demonstrated. It can induce a sizable change of 
the cross section through the significant dependence of the stau masses and 
maxing as well as the h°T^r^ couplings on the phase. At e + e~ colliders with 
cm. energy around yjs = 1 TeV and with very high luminosity J Cdt ~ 500 
fb -1 , the typical cross section is of the order of 0(0.002 ~ 0.03) fb depending 
mainly on the universal scalar fermion mass mo, in particular for light m 
(m < 150 GeV), the cross section can exceed the level of a 0.01 fb, leading 
to several dozen events in a few years. In the case where the tau slepton de- 
cays into a neutrino and a chargino, the final state topology should be easy 
to be seen experimentally owing to the clean environment of these colliders. 
Therefore, all of these results can be expected to be detectable at future 
e + e~~ linear colliders. Analyzing these final states will allow to measure this 
important coupling and determine the value of the CP violating phase, thus 
opening a window to probe directly some of the soft-SUSY breaking scalar 
potential and also to shed some light on the SUSY breaking mechanism. 
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